Quantum phase transitions of polar molecules in bilayer systems 
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We investigate the quantum phase transitions of bosonic polar molecules in a two-dimensional 
double layer system. We show that an interlayer bound state of dipoles (dimers) can be formed when 
the dipole strength is above a critical value, leading to a zero energy resonance in the interlayer s- 
wave scattering channel. In the positive detuning side of the resonance, the strong repulsive interlayer 
pseudo-potential can drive the system into a maximally entangled state, where the wavefunction is 
a superposition of two states that have all molecules in one of the two layers and none in the other. 
We critically discuss how the zero-energy resonance, dimer states and the maximally entangled state 
can be measured in time-of-flight experiments. 



Introduction: Systems of ultracold atoms have become 
one of the most promising systems to observe strong cor- 
relation effects in many-body physics. Recent progress 
in the trapping and coohng of chromium atoms [1] and 
polar molecules Q further opens new directions for in- 
vestigating quantum many-body states resulting from 
the anisotropic dipole-dipole interaction The long 
range nature of dipole interaction also makes it possible 
to study physics in spatially separated multi-component 
systems, which have been extensively studied in several 
important subfields of solid state physics: for example, 
condensation of excitons in bilayer quantum well system 
P|, interlayer fcrromagnctism in bilayer quantum Hall 
systems [5j and Coulomb drag in coupled quantum wires 
[6| etc. Therefore it is interesting to investigate what new 
physics one may expect in the similar systems of polar 
molecules. Recent example is the proposed chaining phe- 
nomena for molecules in a stack of 2D traps which 
resembles particle aggregation in colloidal fluids Q . 

In this paper we investigate the quantum phase tran- 
sitions of cold polar molecules trapped in a 2D double 
well potential 0]. The electric dipole moment (D) is 
aligned perpendicular to the layer {x — y) plane by a DC 
electric field (see Fig. (Ha)) so that the system prop- 
erties is controlled by a dimensionless dipole strength, 
Uq = mD^/h^d, with m being the molecule mass and 
d being the layer separation. We find three phases that 
can be observed in three different regimes of Uq: For 
weak dipole strength, Uq <^1, the ground state is just a 
coupled superfluid (Fig. [Ha)). When Uq is increased to 
be above a critical value, Uq ^ 0.71, molecules in differ- 
ent layers can form interlayer bound states, driving the 
system to be a superfluid of dimers (Fig. [I{b)). Finally, 
if the molecules are cooled in the large Uq regime and 
the dipole moment is reduced toward Uq adiabatically 
from above, we demonstrate that the repulsive interlayer 
pseudo-potential can drive the system to a maximally en- 
tangled state as Uq Uq* ~ 1.4, breaking a global C/(l) 
symmetry via a second order transition. Such maximally 
entangled (ME) state is a superposition of two macro- 
scopical states (or called GHZ state [10]) that have all 
molecules in one layer and none in the other (Fig. [ijc)) 
[m , and therefore will not have any interference pattern 
even in a single shot time-of-flight measurement. 
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FIG. 1: Three many-body states we consider in this paper: 
(a) coupled superfluid state, (b) superfluid of dimers, and (c) 
maximally entangled state. 



Pseudo-potential: We start from the low energy scat- 
tering properties between two molecules via dipole inter- 
action. In the strong confinement regime, we can first 
assume only the lowest subband of each layer is occupied 
and no single particle tunneling between them. The 2D 
Schrodinger equation in the relative coordinate can then 
be written to be 

- ^ + 9y) m + ^o/i(r)^(r) i?^(r) (1) 

where Vo/i(r) = / dzidz2 \^Po{zi - d/2)|^ Wo{z2 T d/2)\^ 
Vd{v,zi — Z2) is the bare interaction for the two 
molecules in the same/different layers. Vd(r, z) = 
£)2 ^2 _ 2^2^ _|_ 2;2)5/2 is dipolc interaction with 
r being the relative coordinate in the x — y plane, (poiz) 
is the lowest confined wavefunction and can be approxi- 
mated by a Gaussian wavefunction of width W {W <C d). 
When finite interlayer tunneling (t) is considered, one 
has to diagonalize the full two-particle two-layer Hamil- 
tonian. Its effect to Eq. (1) can be shown to be the 
order of t'^ /{D'^ /W^), and hence negligible in the strong 
confinement regime as we considered here. 

Now, using the standard scattering theory [l2|, we can 
derive the following 2D pseudo-potential: 
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tan4°^/^^^(fc) -5(0, (2) 



which reproduces the same s-wave phase shift 
^j(o)/(i)^^^-j as the bare interaction, V^/i(r), in large 
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FIG. 2: Intralayer (dashed lines) and interlayer (solid lines) 
scattering phase shift in the s-wave channel as a function of 
momentum kd near resonance ((7o = 0.7). The layer width, 
W = O.ld, is set much smaller than the interlayer distance d. 
Note that phase shift (5o) of the interlayer scattering is posi- 
tive, and its sign is reversed for the convenience of comparison 
with the intralayer results. Inset: Results for t/o = 0.1. 



distance, k is the magnitude of the incoming relative 
momentum. We can show that contributions from higher 
angular momentum channels can be safely neglected 
since the typical length scale of dipole interaction, 
mD^ /h^{^ 1.5/im for D ~ 1 Debye and m ^ 100 a.m.u.) 
is much smaller than the typical condensate size. 

Zero energy resonance: In Fig. [21 we show the calcu- 
lated s-wave phase shift as a function of kd for C/q — 0.7 
and Uo = 0.1 (inset). When the dipole strength is 
weak(inset), the phase shift of the interlayer scattering 
is always much smaller than that of intralayer one as ex- 
pected, but it becomes much larger when Uo is larger. 
In Fig. [21 we show the numerically calculated low en- 
ergy {kd 0) s-wave phase shift as a function of Uq, 
and find a resonance at Uq = Uq ^ 0.71. Similarly 
to the Feshbach resonance in typical cold atom systems 
, this zero energy resonance is due to the formation of 
an interlayer bound state (dimer). Interaction between 
two dimers can be obtained by integrating out the dimer 
wavefunction (not shown here). The system ground state 
for Uq > C/q is a superfluid of dimers with a finite bind- 
ing energy (about 0.1 x h'^/md'^ at Uq — 1), but the 
actual phase transition position might be shifted from 



Uq due to the interaction between dimers. We also note 
that the quantum phase transition from a coupled super- 
fluid to dimer superfluid near Uq is belong to the Ising 
type transition, because the dimer superfluid phase just 
breaks a U{l)/Z2 symmetry, similar to bosnoic systems 
near Feshbach resonance as discussed in Ref. [lj|. We 
point that these result cannot be reproduced even quali- 
tatively within the Born approximation in the literature 
, which is valid only when the dipole strength is very 
weak {Uq < 1). 

Condensate size near resonance: It is interesting to 
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FIG. 3: Zero-energy resonance of the interlayer (solid line) 
scattering phase shift as a function of Uo- Dashed line is 
for the iniralayer scattering. Inset: Calculated condensate 
radius as a function of Uo- Solid and dashed lines are for 
harmonic oscillator length, aho = A/ft/moJn = lOd, and 5d 
respectively. Number of molecules in each layer is A'^ = 10^ 
and other parameters are the same as in Fig. (2] Our meanfield 
treatment of Eq. ([3} may fails near resonance and therefore 
we eliminate the data between the two dotted lines. 



study how the condensate size is changed when the dipole 
strength Uq is tuned across the resonan ce p oint. Us- 
ing a Gaussian variational wavefunction [i6|, 5'±(r) = 
e~l'"l^/^^^(/3o(^ Td/2), for the condensate wavefunc- 
tion in the upper(-f-) and the lower(-) layers, the radius 
R in the negative detuning side is then obtained by min- 
imizing the following meanfield energy: 
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Here t^n is the in-plane trapping frequency, and 
4°)/(^'(i?) = -/o°°xtan4°)/(^^(x/i?)e--V4 are the di- 
mensionless interaction energies, t is single particle tun- 
neling amplitude and here contributes to a constant only 
within meanfield approximation. We can also apply a 
similar method to describe the condensate size of dimers 
in the positive detuning side {Uq > Uq). When Uq is 
well above Uq (i.e. large binding energy of dimers), 
the low energy scattering does not break a dimer and 
the phase sift can be calculated from the interaction be- 
tween dimers after taking into account their bound state 
wavefunction [13. In the inset of Fig. [3 we show the 
calculated condensate radius as a function of Uq. In 
the negative detuning side, the condensate size decreases 
gradually as Uq approaching Uq from below due to the 
increasing attractive interlayer pseudo-potential (main 
plot). On the other (positive detuning) side, the size 
of the dimer condensate grows rapidly due to the repul- 
sive interaction between dimers. Although the meanfield 
calculation may not be reliable when very close to the 
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resonance regime due to the strong momentum depen- 
dence of interlayer pseudo-potential, it is reasonable to 
expect that the sharp shrinking of condensate size near 
Uq is still qualitatively true. Therefore measuring the 
dramatical change of condensate size near resonance can 
provide a clear evidence of zero energy resonance as well 
as the dimer state in the bilayer system. 

Maximally entangled state: In the previous discus- 
sion we concentrated on the situation where the dipole 
strength Uq is initially small and adiabatically increased 
to be above Uq. However, in a realistic experiment, the 
electric dipole moment can be so strong that molecules 
are cooled directly in the large Uq regime with very small 
transition rate to the dimer state. It is therefore in- 
teresting to study how the many-body metastable state 
is changed when the dipole strength Uq is adiabatically 
tuned toward the critical value (Uq) from above. From 
Fig. [3] one can see that there are two regions of interest 
in this positive detuning side: one is for Uq > Uq* ^ l.A 
where the effective interlayer interaction is repulsive but 
still smaller than the intralayer interaction in the long 
wavelength limit, and the other is for Uq < Uq < Uq* 
where the interlayer pseudo-potential is larger than the 
intralayer one. Since both inter- and intra-layer interac- 
tions are repulsive for Uq > Uq, hereafter we may neglect 
the in-plane trapping potential and consider a homoge- 
neous system for simplicity. Assuming all dipoles are in 
the zero momentum state at zero temperature, we can 
write the following effective Hamiltonian 
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where 60(6^) are boson creation operators in the up- 
per/lower layer at k — with and ni, being their num- 
ber operators, gi e 
(i — 0) or inter- (i 



tan 5q'^ {k 0) is the intra- 
1) meanfield energy with being 
the condensate area. Ag = gi ~ go < for Uo > Uq* and 
Ag > for [/* < Uq < U^* . In the second hue of Eq. g]) 
we have used fia+nb = 2N as the conserved total number 
of molecules. We note that although Eq. ([4]) looks sim- 
ilar to the two-site Bose-Hubbard model with inter-site 
interaction, the physics described by Eq. (4) is different 
from the well-known superfluid to Mott-insulator phase 
transition [18]. For example, the lowest excitation state 
of our system is always the in-plane gapless phonon mode 
and therefore no charge gap or commensurate filling are 
expected even when t is reduced to zero. 

Before showing the calculation results, it is instruc- 
tive to discuss the analytic solutions in three different 
limits: Firstly, in the limit of t <C jAgj and Ag < 0, 
the ground state wavefunction is very close to a Fock 
state: I'I'Fock) = mOo^^o^IO)) with almost zero inter- 
layer phase correlation (or the phase stiffness is very 
weak). Secondly, for Ag = but with finite tunnel- 
ing, the ground state is a condensed symmetric coher- 
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FIG. 4: Many-body wavefunction. Cm, for total number of 
2N — 200 dipoles. Solid, dashed, dotted, and dash-dotted 
lines are for Ag/Nt = —3, 1.05, 1.2, and 3 respectively. Inset: 
particle number variation as a function of Ag/t. solid, dashed, 
and dotted lines are for dipole number — 40, 100, and 500 
in each layer. 
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condensate are now phase locked by the single particle 
tunneling so that there will be a true phase correla- 
tion, which can be measured in series of time-of-flight 
experiments. Finally, in the limit of Ag 3> i > 0, 
the total energy is minimized by (hanb) = 0, i.e. all 
dipoles are in one of the two layers and none in the 
other. The most general ground state wavefunction is 
a superposition of two macroscopic states, I^'me) = 

^ (cos (al^^ + sin ^e'^bl^^^ |0>, with tilted angle 



^ and phase x being arbitrary. Such state is also known 
as a kind of Greenberg-Horne-Zeilinger state [l3|, which 
maximizes the entanglement in many measures. Note 
that the maximally entangled (ME) state we consider 
here is spontaneously generated as an exact eigenstate of 
the system, stabilized by the many-body effects. 

The Hamiltonian of Eq. ([4]) can be eas- 
ily diagonalized in the Fock states basis: 



{m\{2N - rny.y^^^ al"\ 



t2Af- 



'|0), where 

TO = 0, 1, • • • , 2N is the number of dipoles in the upper 
layer. The ground state wavefunction can therefore be 
written to be \'^g) = C'm|</>m) with {C^} being 

the eigenvector associated with the lowest eigenenergy of 
Hamiltonian matrix, (0m|-ff|0n)- Similar approach can 
also be applied to systems of finite trapping potential. 
In Fig. [4] we show the exact numerical results of the 
ground state wavefunction (Cm) for different values of 
Ag/t. One can see that for smaller Ag/t (solid line), the 
wavefunction is peaked at to = with a finite distribu- 
tion width to gain tunneling energy. When Ag/t is close 
to one, the single peak distribution becomes unstable 
and tends toward a double-peak distribution. Increasing 
Ag further (i.e. reducing Uq in the positive detuning 
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side) drives the distribution to peak near m = and 
27V, indicating an ME state as discussed above. In the 
inset of Fig. [U we show the particle number variation 
of the ground state, {AN"^) = (V'cKaJao " 4^o)^I*g), 
as a function of Ag/t for various numbers of dipoles per 
layer {N). One can see that in the thermodynamic limit 
(i.e. keeping Ag cx iV/fi fixed as ^ oo), there is a 
sharp phase transition exactly at At; /t ~ I, above which 

\J {AN"^) /2N becomes finite and eventually saturates. 

To understand such a sharp phase transition from a co- 
herent state to the ME state, we can rewrite Eq. ([4]) as a 
spin model (upto a constant) [ll[: H = —2tSx — ^S^, 
where = 5(0^60 + 4"o); = ^{b\af^ - alK)^ 
and Sz = 5(^0*^0 ~ ^0^0) • The total spin is then 
given by = \{fia + ni,){na + fib + 2). In the ther- 
modynamic limit [N 00), we can treat S a clas- 
sical spin: S = iV(sin6' cos (/), sin6' sin </), cos 0) with 9 
and (p being the spin angles in 3D space. Therefore 
the ground state is obtained by minimizing the energy 
E{9,(j))/N — ~2tsm9 coscj) — Agcos'^9 with respect to 
9 and (f). Since must be zero to gain the tunneling 
energy, we can expand E{9, 0) to the leading order of 
77 = — 7r/2 and obtain a Ginsberg-Landau type energy: 
E{r],0)/N = -2t + (t- Ag)r]^ + ^{Ag - i/4)V: which 
shows a clear second order phase transition at Ag = t. 
The variation of particle number (i.e. variation of Sz) 
scales as {Ag — 1)^^^ near the transition point. 

Before concluding, we remark on several experimental 
issues for observing the maximally entangled state in the 
bilayer system. First, for a typical polar molecule Uq can 



be as large as 4-5 and can be easily reduced to zero by 
decreasing the external DC electric field. Secondly we 
can show that phase separation (i.e. dipoles accumulate 
inhomogeneously in layers) is unlikely to occur because 
it causes extra kinetic energy compare to the homoge- 
neously entangled state. Thirdly, the three-body colli- 
sions induced transition to dimer states can be strongly 
suppressed in the ME state, because most molecules now 
are in one of the layers and very few molecules are in the 
other layer. As a result the ME state we proposed here 
should be a long-lived meta-stable state and hence can 
be easily observed in experiments. Finally, unlike the in- 
terference pattern of two independent condensates p^ . 
the fringe contrast of such entangled state will disappear 
even in a single-shot time-of-flight measurement as Uo 
is adiabatically tuned to be lower than Uq* from above. 
This is because the maximally entangled state, I'^mf), 
is a superposition of two macroscopic states, and hence 
very fragile to collapse in any quantum measurement. 
Therefore the disappearance of interference pattern in 
the positive detuning side could be a direct experimental 
evidence of such maximally entangled state. 

In summary, we demonstrate that loading polar 
molecules into a bilayer system can result in several in- 
teresting new physics, including zero energy resonance, 
interlayer bound states (dimers), and a second order 
quantum phase transition toward a maximally entangled 
state. These new phenomena should be easily observable 
using present experimental techniques. 
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